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8 نشریه پژوهشهای علوم و صنایع غذایی ایران 


6۲ صصا. زتا و1۴ //:۵9)) 


مقاله پژو هش 


جلت ۲ شماره ۶۲ قرو فیی‌سا خسیشفت. ۱۴۰۳ من ۸۵-۹4 


مدل‌سازی انتقال حرارت در کنسرو حاوی عصاره مالت در سته‌بندی نیمه سخت بر پایه 


الومینیوم 


وه 7 ۲ ۰ ۳ 2 ۳ 
سعیده حمیدی - نفیسه زمین‌دار 8 - نیره قلی پور شهر کی 
تاریخ دریافت: ۱ 


تاریخ پذیرش: ۱۴۰۱/۱۲/۱۵ 


چکیده 

در این پژوهش مدل عددی انتقال حرارت در کنسرو عصاره مالت با بریکس ۶۰ در بسته‌بندی آلومینیومی نیمه‌سخت با نرم‌افزار فلوئنت توسعه داده شد. شکل 
هندسی کنسرو حاوی نمونه» توسط نرم‌افزار گمبیت رسم و شبکه‌بندی مناسب با فاصله گره‌های ۰/۱6۲ و ۰/۲۵۲0 اعمال شد. سپس انتقال حرارت درون کنسرو 
در طول فرآیند حرارتی با استفاده از نر‌افزار فلوئنت شبیه‌سازی شد. پروفیل‌های دمایی فرآیند حرارتی. شکل و محل ناحیه کند حرارتی در کنسرو نمونه بررسی شد. 
همچنین تآثیر وجود سرفضا (فضای خالی بالای ظرف) بر انتقال حرارت مورد بررسی قرار گرفت. شبکه‌بندی مناسب برای شبیه‌سازی» شبکه با فاصله گره‌های 
۰ بود. نتایج شبیه‌سازی نشان داد که محل ناحیه کند حرارتی در کنسرو عصاره مالت دارای سرفضا در پایان مرحله حرارت‌دهی در محدوده ۳/۲۷۵۲۲- 
۰۳/۳۷ ۱۷۰/۳۵0۲۵ ۰/۰۲>۲>۰/۸0۲۲- و در مدل فاقد سرفضا در پایان مرحله حرارت‌دهی در محدوده ۳/۴۶۳/۰۵۵ ۰/۴۸6۲۲> 
۴۷ ۵ ۳/۵۸>>۳/۷۶- است. بین دماهای پیش‌بینی شده و دماهای حاصل از تکرارهای آزمایشگاهی در سطح احتمال ۱ تفاوت معنی‌دار مشاهده 
نشد. پروفیل‌های دمایی حاصل از شبیه‌سازی در حالت‌های دارا و فاقد سرفضا در مراحل مختلف فرآیند حرارتی دارای مشابهت قابل قبولی بودند. 


واژه‌های کلیدی: بسته‌بندی آلومینیومی نیمه انعطاف‌پذیر دینامیک سیالات محاسباتی» عصاره مالت» ناحیه کند حرارتی 


مقدمه بردن اسپورهای بیماری‌زا در محصول غذایی است. در صنعت کنسرو 
اغلب عملیات‌های فرآوری مواد غذایی جهت افزایش عمر مفید این فرایند معمولا با استفاده از بخار برای گرم کردن غذا تا یک دمای 
محصول از طریق حذف فعالیت نامطلوب میکروبی طراحی شدهاند. 


بیش از یک قرن است که فرآیند حرارتی در تولید غذاهای کنسرو شده 


مشخص و سپس نگهداری غذا در آن دما برای یک بازه زمانی کافی 
به‌منظور نابودی میکروارگانیسم‌ها انجام می‌شود. در استرلیزاسیون 


با عمر ماندگاری بالا نقش بسزایی دارد. شناخته شده‌ترین عملیات 
فرآوری که در آن انرژی حرارتی به کار می‌روده استرلیزاسیون است. 
این فرآیند شامل استفاده از حرارت در دما و زمان خاصء» جهت از بین 


حرارتی» معمولا دو روش اعمال می‌شوده ناپلشت کردن" که در آن 
محصول غذایی قبل از بسته‌بندی استریل می‌شود و روش دیگر قوطی 


(خوراسگان» اصفهان, ایران 
(#- نویسنده مسئول: 11151.26.1۳ ۵) عل صتصح 1.2 :اتقحصط) 


۳- دانشجوی دکتری علوم و صنایع غذایی» واحد شهر کرد دانشگاه آزاد اسلامی» شهر کرد ایران 
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استریل می‌شود ( ,.اه ۶ اعد :2004 متصعطم اب۸0 ک ۳۵۳۱0 
2 قوطی‌های کنسرو پس از تولید باید برای مدت‌های طولانی 
قابل نگهداری باشند. برای این منظور دما و زمان سالم‌سازی باید تا 
حدی باشد که مقاوم‌ترین میکروارگانیسم‌های احتمالی موجود در 
محتوای قوطی‌ها نابود شوند. در کنسروسازی زمان و دمای 
استریلیزاسیون بر اساس مقاومت دمایی اسپور کلستریدیوم بوتولینوم! 
محاسبه می‌شود. اما به‌طور کلی زمان» درجه حرارت پخت و 
استریلیزاسیون در مورد هر محصول و هر نوع بسته لازم است جداگانه 
مورد محاسبه قرار گیرد (2002 ,۸ ۶ تصعطت 001ه). مدت زمان 
لازم برای استریلیزاسیون محصول به ویژگی‌های محصول, نوع و اندازه 
بسته و نحوه قرارگیری در اتوکلاو و محیط حرارتی بستگی دارد (» 
9 ,۸ 0۵ آبا۸). فرایند استریلیزاسیون علاوه بر پهبود عمر 
انبارمانی غذا و ايمن ساختن آن از طریق تخریب میکروارگانیسم‌های 
بیماری‌زاه می‌تواند بر کیفیت تغذیه‌ای محصول از جمله ویتامین‌های آن 
اثر داشته باشد. بنابراین یکی از چالش‌های فرایند استریلیزاسیون 
کنسروهاء تولید محصولی سالم» ايمن و در عین حال دارای کیفیت 
تغذیه‌ای مطلوب است (2002 ,.۵1 2 تصعطت) ۸001). حرارت اضافی 
بر کیفیت و ارزش تغذیه‌ای محصول اثر منفی می‌گذارد. تخمین سرعت 
انتقال حرارت برای رسیدن به شرایط بهینه فرایند و بهبود کیفیت 
محصول ضرورت دارد. درک و فهم بهتر فرایند حرارتی منجر به بهبود 
انجام فراینده صرفه‌جویی در مصرف و ذخیره انرژی می‌شود ( .41 1» 
8 ,۵ 6 ۴۵02۷1 :1999 تصعطن آبل‌دش). بنابراین باید مطالعاتی 
در راستای مدل‌سازی» پیش‌بینی و تعیین الگوی حرارت‌دهی مواد 
غذایی. در طی فرآیند حرارتی صورت گیرد تا بتوان طراحی بهتری 
برای فر آیند استرلیزاسیون داشت. به این منظور از علم دینامیک سیالات 
محاسباتی " به‌عنوان ابزاری قدرتمند استفاده می‌شود. دینامیک سیالات 
محاسباتی روشی برای شبیه‌سازی. جهت حل مسائل جریان سیال و 
انتقال حرارت است که در آن معادله پیوستگی, معادله ممنتوم و معادله 
بقای انرژی به‌صورت عددی حل می‌شوند تا سرعت دماء برش " فشار 


و دیگر پارامترهای موجود در یک سیستم جریان سیال پیش‌بینی شوند 
( ک ۵۷1 :2021 باه ۶ طفلههند۲ :2006 مصتاگ عک طماع0( 


1 ,016). فرآیند تشکیل مالت شامل جوانه زدن دانه‌های سالم 
جوء تحت شرایط خاص از نظر دماء رطوبت و زمان مشخص می‌باشد. 
به‌نحوی که اندوسپرم دانه (تأمین کننده رشد جوانه است» تحت 
واکنش‌های بیوشیمیایی و تجزیه آنزیمی قرار می‌گیرد. در اين حالت 
ویتامین‌های گروه ب و آنزیم‌ها بویژه آنزیم‌های تجزیه کننده نشاسته 


۱ 0051۳۱1۱۱ -1 
(ناتت) ومتصمط نا متا آلمطمتادانام‌ططمن -2 


ه ضا ی الق ارات ایتک ی زا سای 
آبی فالههای غلات. خصوضا جو کذ ق آیند. تشکیل مالت: را گذرانده 
بدست می‌آید (6960 ۱۷/312 051121. این محصول به‌طور 
گسترده در صنایع غذایی و همچنین دارویی مورد استفاده قرار می‌گیرد. 
استفاده عمده از مالت تولیدی در جهان در صنعت آبجوسازی است. با 
توجه به اينکه کار با بسته‌های آلومینیومی نیمه‌سخت راحت است و از 
آنجایی که بسته‌های خالی از این نوع به راحتی در یکدیگر فرو می‌روند. 
فضای کمی را اشغال می‌کنند و بنابراین حمل و نقل آن‌ها آسان است 
لذا در این تحقیق از بسته‌های نیمه‌سخت با پایه آلومینیوم برای کنسرو 
کردن استفاده شد. با توجه به اینکه بسته‌های آلومینیومی نیمه‌سخت 
اخیراً بهبازار مصرف معرفی شده‌انده اطلاعات کمی در مورد نحوه توزیع 
دما در طول استریلیزاسیون آن‌ها وجود دارد. در اين پژوهش با توجه به 
ویژگی‌های ترموفیزیکی اندازه‌گیری شده و با استفاده از نرم‌افزار فلوئنت 
مدلی ریاضی برای توصیف انتقال حرارت در بسته‌های آلومینیومی نیمه 
سخت در طی فرایند حرارتی عصاره مالت به‌عنوان تابعی از دمای 
اتوکلاو توسعه داده شد و پروفیل‌های دماء شدت فرآیند حرارتی و 
کندترین ناحیه حرارتی " قوطی پیش‌بینی شد. سپس داده‌های پیش‌بینی 
شده به‌وسیله مدل با مقادیر آزمایشگاهی مقایسه شد تا اعتبار مذل 


بررسی گردد. 


مواد و روش‌ها 
تهیه نمونه 

در این پژوهش عصاره مالت با درجه بریکس ۶۰ تهیه شده در 
کارخانه به مالت شهرکرده جهت مدل‌سازی فرآیند انتقال حرارت مورد 
استفاده گرفت. 


تجزیه شیمیایی نمونه 

تجزیه شیمیایی عصاره مالت شامل اندزه‌گیری پروتئین به روش 
کلدال, چربی به روش ورنره محتوای خاکستر و رطویت ( 151181 
7 ۱1/11312). کربوهیدرات (1952 ,10005 و محتوای فیبر 
خام در سه تکرار انحام شد. 


تعیین خواص ترموفیزیکی 
چگالی 

جرم حجمی عصاره مالت در محدوده دمایی )۲۵ تا ۸۰۳6 در سه 
تکراراندزه‌گیری شد. جهت اندزه‌گیری چگالی نمونه مورد نظر به‌دلیل 


۲ -3 
عم ۲6۵1182 10۷76۵۲ -4 


حمیدی و همکاران» مدل‌سازی انتقال حرارت در کنسرو حاوی عصاره مالت در بسته‌بندی نیمه‌سخت بر پایه آلومینیوم ۸٩‏ 


لزجت بالا و محبوس شدن هوا در آن. نمونه‌ها به‌مدت ۲ ساعت تحت 
خلاً قرار گرفتند (2015 .۵1 6۶ طحطلصها ۷). 


لزجت 

ندزهگیری لزجت نمونه توسط ویسکومتر بروکفیلد مدل 0۷-111 
5۴ ساخت آمریکا در بازه‌های ۱۰ درجه‌ای از ۳۰۳6 تا ۸۰۳6 انحام 
شد (2012 ,نت۸۱ > 6222۷71 


گرمای ویژه 
گرمای ویژه نمونه با استفاده از روابط تجربی ارائه شده توسط چوی 
و اوکوس و معادله (۱) که بر پایه درصد اجزای تشکیل دهنده نمونه 
است در محدوده دمایی ۲۰۳۵ تا ۸۰۳ محاسبه شد ( عک طنطه6 
6 را 0), 
)۱ نون دحایظ < وب) 
در این معادله 0۳1 گرمای ویژه جزء سازنده 1 در فشار ثابت و :> 
جزء جرمی ماده سازنده 1 است. 


ضریب هدایت حرارتی 
ضریب هدایت حرارتی نمونه با استفاده از روابط تجربی ارائه شده 
توسط چوی و او کوس 9 معادلات ۲ 9 ۳ در بازه دمایی تا ۸۰۲ 


0/7 ۷ 
۷ ( 0 ِِ« 
)۳ + 5 وک 


ذام و رم چگالی جزء 1 ام ((/18) است. 


حل مسئله با استفاده از نرم‌افزار فلوئنت 


۱ 
«۹ 


برای بررسی عددی و توصیف انتقال حرارت در کنسرو» ابتدا هندسه 
ظرف توسط نرم‌افزار گمبیت طراحی و توسط همین نرم‌افزار شبکه‌بندی 
شد. پس از فراخوانی فایل شبکه‌بندی شده در نرم‌افزار (فلوئنت) حل 
مسئله انجام شد. 


طراحی هندسه مدل 

در این مطالعه دو مدل هندسی دارا و فاقد سرفضا طراحی شد. لازم 
به ذکر است منظور از مدل فاقد سرفضا این است که در طراحی مدل 
سازی انتقال حرارت در نرم‌افزانه فقط بخش ماده غذایی در نظرگرفته 
شد و بخش سرفضا برای مدلسازی لحاظ نشد و شرایط مرزی در بخش 
فوقانی ماده غذایی مشابه شرایط دیواره جانبی درنظر گرفته شد. 

در طراحی مدل ابتدا یک مربع به ضلع ۸/۶۰6۰ رسم شد و در 
فواصل ۲۲۰ از گوشه‌های آن کمان‌هایی در نظر گرفته شد. سپس 
قسمت‌های اضافی حذف شده و خطوط موجود به صفحه تبدیل شد. 
سپس صفحه ایجاد شده در راستای منفی محور 2 به اندازه ارتفاع ظرف 
(۴/۱۵۵۲) و با زاویه ۵ درجه گسترش داده شد. در طراحی مدل دارای 
سرفضا مراحل بیان شده جهت طراحی سرفضا با ارتفاع ۰/۴۰۲۲ انجام 
شد و سپس شکل ایجاد شده در راستای منفی محور 2 به اندازه ارتفاع 
ظرف (۳/۷۵۵) و با زاویه ۵ درجه گسترش داده شد تا شکل هندسی 
ظرف مطابق شکل ۱ بدست آید. 


شبکه‌بندی 

فضای محاسباتی باید ثبکه‌بندی شود تا نرم‌افزار تحلیل گر بتواند 
برای واحدهای حاصل» حل عددی انحام دهد. به‌منظور شبکه‌بندی 
مدل, ابتدا شبکه‌بندی روی سطوح و سپس شبکه‌بندی روی حجم 
صورت گرفت. در طراحی مدل یک بار از شیکه‌بندی عته ۰/۱ و بر 
دیگر جح ۰/۲ استفاده شد. شکل ۲ مدل شبکه‌بندی شده با اندازه شبکه 
۲ را نشان می‌دهد. 


شکل ۱- شکل هندسی ظرف (سمت راست: دارای سرفضاء سمت چپ: فاقد سرفضا) 


ده ۵20 )م۱1 :میا دص 0۵20 )۱ :)2 همع مصنداری صسصتصساه مزوتدنصعی :1 .۲12 


۰ . نشریه پژوهشهای علوم و صنایع غذایی ایران. جلد ۰۲۰ شماره ۱. فروردین-اردیبهشت ۱۴۰۳ 


شکل ۲- مدل شبکه بندی شده با اندازه شبکه ‏ ۰/۲ 
(صه 0.2) هجاوم 6 ۵۶ ۵۵061 31 من مم‌طاومصه 160/۳۲۲۵۲۱۵ 2۰ ۲۱۵۰ 


تعیین مرزها در نرم‌افزار گمبیت 

در این مرحله کلیه مرزها و نوع مواد از نظر سیال یا جامد بودن در 
نرم‌افزار گمبیت تعریف شد و در انتها مدل طراحی شده بصورت فایل 
شبکه‌بندی شده به نرم‌افزار تحلیل‌گر (فلوئنت) ارسال شد. مسئله مورد 
نظر به‌صورت حل با دقت دو برابر و در حالت سه بعدی حل شد. 


خواندن و بررسی شبکه 

مدل طراحی شده نهایی پس از فراخوانی در نرم‌افزار تحلیل گر 
(فلوتنت) ابتدا مورد بررسی کیفی" قرار گرفت. با توجه به اینکه هندسه 
مسئله در مقیاس سانتی‌متر تعریف شده است» تمامی مقیاس‌های مورد 
استفاده بر حسب واحد سانتی‌متر در نظر گرفته شد. 


انتخاب شیوه محاسباتی و فرمول‌بندی حل مسئله و تعریف 
معادلات حاکم 

در حل مسئله از شیوه پیوسته ضمنی استفاده شد. معادلات حاکم 
با شرایط فشار پایه » به‌صورت ضمنی" در حالت گذرا" و معادلات 
غمتن در مویخلی تفه عدل فیزیکی مهرد نظر بر اه جریا آز 
پایه‌گذاری شده است. به‌منظور بررسی انتقال حرارت معادلات انرژی " 
فعال شدند. شتاب ثقل برای محور ۰۷ 10/52 ۹/۸- در نظر گرفته شد. 


تعیین خواص ترموفیزیکی 


1- 0 

2- ۳۲۵۵۹۱9۲ 0 

ازمنامصص1 -3 
07 -4 

5- 2۴ 

کممتاها۲0 و۲6۲ -6 


تعبین قترابط موزیق اولبه 

شرایط مرزی" با توجه به شرایط واقعی اعمال شده بر دیواره‌های 
ظرف در نرم‌افزار تعریف شد. به این صورت که برای دیواره‌های جانبی» 
درب و کف ظرف ویژگی‌های ممنتوم به صورت دیواره ثابت" و ویژگی - 
های برش به‌صورت بدون لغزش* در نظر گرفته شد. شرایط مرزی 
حرارتی مربوط به مدل فاقد سرفضا؛ برای دیواره‌های جانبی» کف و 
درب ظرف به صورت دمای ثابت و منطبق با پروفیل حرارتی محیط 
داخلی اتوکلاو (مطابق با شکل ۳) تعریف گردید. شرایط مرزی حرارتی 
مربوط به مدل دارای سرفضا؛ برای دیواره‌های جانبی و کف ظرف مشابه 
مدل فاقد سرفضا در نظر گرفته شد. اما با توجه به در نظر گرفتن سرفضا 
در این مدل؛ شرایط مرزی حرارتی برای قسمت فوقانی (درب) ظرف به 
صورت انتقال حرارت از نوع همرفت و براساس ضریب همرفت بخار 
تعریف گردید. خصوصیات مربوط به ظرف نیز براساس جنس آلومینیوم 
و با ضخامت دیواره 0 ۰/۰۴ تعریف شد. شرایط اولیه همان دمای 
اولیه محصول در نظر گرفته شد. 

جهت دستیابی به همگرایی بهتر؛ فاکتورهای مادون رهایی " 
کوچکتر یا مساوی یک در نظر گرفته شد. همچنین فشار؛ استاندارده 
معادلات انرژی و ممنتوم از نوع درجه دو بالا دست" و الگوریتم مربوط 
بگوین سخ کشا سس سس اتف انم فیرمدات 
ظرف به‌عنوان نقطه شروع محاسبات تعریف و در نهایت حل مسئله راه 
تا نت سقایر یرای بای تاقماندفهای وگن رعاش 
راستای محورهای >2 و۷ و 2 مقدار ۱۰۳ و برای انرژی ۱۰ در نظر 
گرفته شد (2015 ,۵1 6 طقط‌لمهاه ۷ 


جمتاتصم تدلصتا0ظ -7 
۵ 210100277 -6 

9- ۱۲0 0 

10- ۲/806۲ 0 

فصو تع0ت0 ٩60080‏ -11 
0 -12 


حمیدی و همکاران. مدل‌سازی انتقال حرارت در کنسرو حاوی عصاره مالت در بسته‌بندی نیمه‌سخت بر پایه آلومینیوم ٩۱‏ 


زمان 


(صتجص) عحطرزً ]" 


دما 


۰ 

۰ ۶۰ 
۰ 
هِ 


همع ]1 


شکل ۳- تغییرات دمای اتوکلاو در طول زمان در فرآیند حرارتی عصاره مالت 


/۱۵۱ 0 ماوعع 0۵ امصمط م6 ما ,و۲ مرها 0۴۲۵۵۲ ومع۵ط .3 ۲۱۵۰ 


انجام فرآیند حرارتی به‌منظور معتبرسازی مدل 

در این مرحله از تحقیق, معتبرسازی پیش‌بینی‌های تئوری توزیع 
دما طی فرایند حرارتی مورد آزمایش قرار گرفت. ظرف مورد نظر با وزن 
استاندارد از نمونه پر و توسط ماشین آلکان" با المنت حرارتی ۲۸۰۳۵ 
درب‌بندی شدء ترموکوپل برای اندازه‌گیری توزیع دما در مرکز هندسی 
ظرف مورد نظر قرار گرفت. ترموکوپل بوسیله کابل‌های ۳1۱۰۰ به 
دیتالاگر " دستگاه ایلب" متصل و ظرف در داخل اتوکلاو قرار گرفت. 
اطلاعات دما و زمان در فواصل یک دقیقه‌ای از طریق نرم‌افزار ۳۲-۷۵1 
ثبت شد. کلیه مراحل معتبرسازی در کارخانه صنایع غذایی چیکا سپاهان 
واقع در اصفهان شهرک صنعتی مبا رکه انجام شد. 


معتبرسازی مدل و ارزیابی خطا 

معتبرسازی مدل براساس مقایسه میانگین داده‌های حاصل از 
تکرارهای آزمایشگاهی (از طریق نصب ترموکوپل در داخل ظرف و ثبت 
اطلاعات دما در طول فرآیند حرارتی به‌دست آمدند) با داده‌های پیش- 
بینی شده توسط نرم‌افزار از طریق آزمون ) انجام گرفت. میزان ريشه 
مربعات میانگین خطا" نیز به‌عنوان معیاری برای برازش از طریق معادله 
(۴) محاسبه شد (2012 ,ته]ت۲ ک با0ع۳۳0۵0). 


1 
(۴) در ههنا هنزو[ ری یت 1۷ :]1 


مصتطمعمه صععآ۸ -1 
6۲ 12/۵ -2 


پایین بودن مقدار ریشه مربعات میانگین خطا نشانگر اعتبار بالای 
مدل است (2007 ,8عطهط۱۷]0 :2004 اه 6 ط0وموصنک). 


محاسبه ارزش استریل‌کنندگی 

در مورد غذاهای کنسروی کمترین 9۳1 ممکن جهت رشد 
میکروارگانیسم کلستریدیوم بوتولینیوم ۴/۶ است. اسپور این 
میکروارگانیسم دارای بالاترین میزان مقاومت گرمایی در میان 
میکروب‌های غذایی بیماری‌زا است. برای انجام فرآیند حرارتی 
استریلیزاسیون جهت غذاهای کم اسید (دارای 11 بالای ۴/۶) فرآیند 
حرارتی می‌باید به گونه‌ای باشد که احتمال حضور اسپور کلستریدیوم 
بوتولینیوم را به ۱۰۳۳ در هر بسته برساند. برای محاسبه میزان ارزش 
استریل کنندگی از معادله (۵) استفاده شده است. 

۵ ۸ و < از 

میزان دمای مرجع (:1) با توجه به دمای اتوکلاو يا دمای پیش 
پیش که عبط سل برایر ۱۲۱/۱۳۵ و رای بان هافر نظار گرفند 
شد (2001 ,۵1 61 نداهعنل. 


اتالیز شیسایی خر کنبات کشکین. دهنله عصاره مات در خقول ٩‏ 


ارائه شده است. 


3- 120 
4- 1۵0] 116286 ۹0۵۲6 ۴۲ 


۹ 


نشریه پژوهشهای علوم و صنایع غذایی ابران. جلد ۰۲۰ شماره ۱. فروردین-اردببهشت ۱۴۰۳ 


جدول ۱- آنالیز شیمیایی عصاره مالت 


۲ ۵1 0۴ وزوباحصه امءتصمطن -1 ما9و 1 


آب کربوهیدرات . پروتئین چربی خاکستر ‏ فیبرخام 
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خواص ترموفیزیکی 

به‌منظور افزايش تأیید مدل پیش‌بینی شده خواص ترموفیزیکی 
نمونه به‌عنوان تابع دما اندازه‌گیری شد. ضریب انبساط حجمی (6) ۴" 
۰ برآورد شد. معادلات (4۶ (۷ (۸) و )٩(‏ به‌ترتیب تغییرات 
چگالی. گرمای ویژه, هدایت حرارتی و لزجت را به‌عنوان توابعی از دما 


نشان می‌دهند. 
۶ 6 + 0.10687 + *0.005971- < ۵0 
)۳ 6 + 9.5757 < م) 
)۸( 7 + 0.00097 < 1 
٩)‏ ۱ 


بررسی روند همگرابی در مدل انتقال حرارت در کنسرو 
عصاره مالت 

بررسی تغییرات دما در بازه‌های زمانی متوالی بر نقطه واقع بر محور 
مختصات به ارتفاع ۲/۷۶ از سر ظرف نشان داد که دما در نقطه مذکور 
به‌صورت تابعی از دمای اتوکلاو ابتدا افزايش سپس کاهش می‌پابد. به 
منظور کنترل روند همگرایی مستله باقیمانده‌های انرژی۱ شکل ۳ و 
ها ار نها ومای ال نر تقط ما کین بورمی قفق که 
هر دو مه تأییدی بر همگرا شدن مسئله هم از نظر باقیمانده‌های انرژی 
و هم از نظر روال منطقی تغیبرات دما بود. 


بررسی اثر شبکه‌بندی متفاوت بر مدل انتقال حرارت در 
کنسرو عصاره مالت 

در مدل طراحی شده از دو شبکه‌بندی ۰/۱۵۳ و ۰/۲۵۲۲ جهت 
حل مستله استفاده شد. دمای محل هندسی قرارگیری ترموکویل په 
مختصات (۲/۷۶ و۰ وه) تا پایان مرحله نگهداری" توسط مدل در دو 
شبکه‌بندی ۰/۱۵ و ه ۰/۲ شبیه‌سازی شد. و توسط آزمون آماری 
مورد مقایسه قرار گرفت. نتایج بدست آمده عدم وجود اختلاف معنی- 
دار در سطح احتمال ۱ را نشان داد. با توجه به نتیجه بدست آمده به 


10و16 زوم -1 
۲0۱ -2 


منظور ساده‌سازی روند حل و جلوگیری از تلاف وقت از شبکه‌بندی 
۶ در ادامه روند تحقیق استفاده شد. 


بررسی پروفایل‌های حرارتی در مدل عصاره مالت در مقاطع 
زمانی مختلف 

پروفیل حرارتی در سطح ۷۰ و ۰7 ۲/۷۶- <2 (از سر ظرف) به 
صورت همدما در طول فرایند حرارتی (حرارت‌دهی و خنک کردن) در 
مدل فاقد سرفضا در شکل ۵ نشان داده شده است. همانطور که در 
شکل قابل مشاهده است دما در داخل ظرف کنسرو بصورت تابعی از 
دمای اتوکلاو ابتدا افزايش و سپس کاهش می‌یابد. همچنین بررسی 
شکل ۵ نشان می‌دهد در مراحل اولیه فرآیند» انتقال حرارت بیشتر به 
صورت هدایت بوده و فرآیند جایجایی نقش کمتری دارد (الف و ب) و 
سپس وارد مرحله گذار از هدایت به جابجایی شده (د و ج) و در ادامه 
فرآیند مکانیسم غالب در انتقال حرارت جابجایی است (ه و و). 


بررسی اثر سر فضای ظرف کنسرو در مدل‌های پیش‌بینی 
شده 

با توجه به این نکته که در شرایط آزمایشگاهی تمامی ظروف پر 
شده از عصاره مالت دارای سرفضا" بودند تأثیر سرفضا بر سرعت انتقال 
حرارت بررسی شد. در این مطالعه ۰/۴۵۳ سرفضا در بخش فوقانی 
تمامی ظروف در نظر گرفته شد. از آزمون آماری ] جهت مقایسه دمای 
محل قرارگیری ترموکوپل در ظرف (۲/۷۶- و۰ و۰) برای مدل دارای 
سرفضا نسبت به مدل فاقد سرفضا استفاده شد. نتایج نشان‌دهنده عدم 
وجود اختلاف معنی‌دار در سطح احتمال ۱ بين دو مدل بود. شکل ۶ 
پروفیل دما در سطح ۷-۰۵۳ را در مقاطع زمانی مختلف در دو مدل 
دارا و فاقد سرفضا نشان می‌دهد. همچنین از آزمون آماری ] جهت 
بررسی تأثیر شبکه‌بندی بر دمای نقطه انتخابی استفاده شد, که نتایج 
حاکی از عدم وجود اختلاف معنی‌دار بین شبکه‌بندی ۰/۱۵ و ۰/۲6۲۲ 


در سطح احتمال ۱/ بود. 


3- ۳۱620/ 6 


حمیدی و همکاران» مدل‌سازی انتقال حرارت در کنسرو حاوی عصاره مالت در بسته‌بندی نیمه‌سخت بر پایه آلومینیوم 


بررسی ناحیه کند حرارتی در مدل عصاره مالت 

در انتقال حرارت به‌صورت هدایت ناحیه کند حرارتی که کمترین 
میزان حرارت را دریافت می‌کند. غالبا در مرکز هندسی توده غذایی قرار 
دارد. ولی مطالعات زیادی وجود دارد که نشان می‌دهد در غذاهایی که 
فقط از طریق همرفت گرم می‌شوندهناحیه کند حرارتی در مرکز هندسی 
قوطی قرار نمی‌گیرد. بلکه در بخشی نزدیک به کف قوطی قرار می 
گیرد. این مسئله هم به‌صورت تجربی و هم با مدل‌سازی کامپیوتری 
اثبات شده است ( ۸01 ی نته۳ :1999 ,اه ۶۶ تصحطن آنط۸ 
4 ,0211)). به‌صورت تجربی و تئوری مشخص شده است که در 
غذاهایی که حرارت در آن‌ها تنها به‌صورت همرفت منتقل می‌شود. 


زارزرززا۸۸۵ ۱ 


تیا 


بررالال/۱۱ 


٩۳ 


محل ناحیه کند حرارتی تغییر می‌کند (2002 ,./۵ 6۶ نصهط اه‌طه). 


همچنین با پیشرفت زمان حرارت‌دهی دمای یکنواخت‌تری در بسته به 
وجود می‌آید و نیروهای جابه‌جا کننده سیال کم می‌شود. در نتیجه 
سرعت جابجایی کم شده و در پایان حرارت دهی جریان چرخشی 
متوقف می‌شود (1999 ,اه ۶ تصعطن ۵۸0۵1 برای بررسی ناحیه 
کند حرارتی و مشخص کردن محل آن؛ برش‌هایی از سطح ۷-۰0۲ در 
مقاطع زمانی مختلف در دو حالت دارای سر فضا و فاقد سرفضا مورد 
بررسی قرار گرفت ([شکل ۶). 


باقی مانده انرژی 
۷ 1۵910121 


شکل ع- باقیمانده انرژی در 00 ۳۲ <] 
صنص 1232 ]2 وافیمنوع۲ ۲۱۵۲۵ .4 ,ع۲1 


خطوط هم‌دمای مربوطه نشانگر این مطلب بود که اگر چه اندکی 
تغییر در ارتفاع و شکل محل قرارگیری ناحیه کند حرارتی در هر دو 
عدل ارام فاقد سرقضا حزمراحل مکسلت فر این گرمانی مفناهده ی < 
شود اما در تمام مراحل فرآیند گرمایی این ناحیه متمایل به یک سوم 
انتهایی ظرف بوده و محل آن ثابت است. بجز در مراحل انتهایی فرآیند 
که این ناحیه متمایل به نیمه بالایی ظرف کنسرو است. در شکل ۵ 
برش‌هایی از سطح ۲/۷۶۵۲۲- 2 نشان داده شد. همانطور که مشاهده 
می‌شود محل و شکل ناحیه کند حرارتی در سطح ۲/۷۶۰- 2 در 
طول فرآیند تغییر می‌کند. به‌طوری که در ابتدای فرآیند حرارتی ناحیه‌ای 
کوچک نزدیک به مرکز مقطع را شامل می‌شود که اين نشانگر غالب 
بودن هدایت حرارتی در اين مرحله از فرآیند است و در اتتهای فرآیند 
حرارتی (مرحله سرد کردن) ناحیه کند حرارتی به دیواره‌های ظرف 
کنسرو متمایل می‌شود. داده‌های حاصل مربوط به تمامی گره‌ها در 
این مرحله هلدینگ» به‌صورت فایل اکسل"» از نرم‌فزار فلوتنت 
اقا ی قانم پیانگر این فکید اسبق که اجه کید حراتی صمل 


1 - ۱۷۲10۲0۹01] ۵ 


طراحی شده کنسرو عصاره مالت دارای سرفضا در پایان مرحله حرارت 
دهی در رک رک 
> >۰/۰۲- و ناحیه کند حرارتی مدل عصاره مالت فاقد 


محدو ده 


سرفضا در پایان مرحله حرارت‌دهی در محدوده ۳/۰۵۵۳0->,/>۳/۴۶- 
۰ ۲۴۴۷۰/۴۸۵ ون ۳/۵۸>۲>۲۳/۷۶- است. داده‌های 
حاصل نشان داد ارتفاع ناحیه کند حرارتی در هر دو مدل دارا و فاقد سر 
فضا دارای مشابهت قابل قبول هستند. 


ار توت کیرات ما وهای مق فا کر فزو گرا 
برای مدل شبیه‌سازی شده و مدل تجربی ارائه شده است. 

همچنین نمودار محاسبات ارزش استریلیزاسیون نسبت به زمان در 
آزمایشگاهی در شکل ۸ نشان داده شده است. 


۴ نشریه پژوهشهای علوم و صنایع غذایی ایران. جلد ۰۲۰ شماره . فروردین-اردیبهشت ۱۴۰۳ 
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شکل ۵- پروفیل دما در مدل فاقد سر فضا در سطح ۷۲۵ ۲- 2 در (الف) من دا (ج) صنص 1۳۲ (ه) صنص ماع 
و در سطح ۰00 ۷ در (ب) صن ۵ (د) هن ۰۱۳۲ (و) صنه ماع 


(0) و(صنجه) 32 < ] 21 (0) زرسنه) 5 <] 26 (2) ۹۳]2»65 2.76- 7.2 10 5۵266 20م )م6 مه نز یمیمص 1 5۰ .۲۱۵ 
(صنصص 60 < ] () رصن 32 < ] (0) زرصتصص) 5 ع ] (ظ) مم‌هسه ۷۶۵ 1 0ج رصتصص) 60 < 1 21 


حمیدی و همکاران. مدل‌سازی انتقال حرارت در کنسرو حاوی عصاره مالت در بسته‌بندی نیمه‌سخت بر پایه آلومینیوم ٩۹۵‏ 
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شکل 71- پروفیل دما در سطح ۷-۰6۳ (الف) مدل بدون سرفضاء (ب) مدل با سرفضا در «۵9-)؛ (ج) مدل بدون سرفضاء (د) مدل با 


سرفضا در « ۲ ۱<؛ (۵) مدل بدون سرفضاء (و) مدل با سرفضا در 


صتصر+زحا 


۳۳۲۸ ؛ (ز) مدل بدون سرفضاء (ح) مدل با سرفضا در 


و() 5 ۲ 24 ۵۵6۵و 6۵0 ۱6 ۱۷۲۵۵ () مرو 0قعظ )بمط6ز۷ ۱۷۱۵06 ره) ممه اه ۷0 صز مسمصی و168۵ 6۰ ۲1۵۰ 
۷۵ () 0266: 6۵0 )۳۷160۵ ۱۷۲۵0۵۵ رع) ورصتصص) 16 < ۲ 2 ۵26و 6۵0 16 ۱۷۱۵0۵۵ (0) عمهوری 0قعظ نامز ۱۷۲۵۵61 0) 
(سنصص) 60 ع ۲ 2۲ 5۵266 620 ۳۲6۵ ۱۷۲۵0۵6 رط) مدمه 0ععظ اسامطز ۱۷۲۵۵6۵ (ع) ز(صتصص) 32 < ۲ 26 5۵266 620 طز ۲۷۲ 


معتبرسازی مدل 

معتبرسازی مدل براساس مقایسه میانگین داده‌های حاصل از 
تکرارهای آزمایشگاهی (از طریق نصب ترموکوپل در داخل ظرف و 
ثبت اطلاعات دما در طول فرآیند حرارتی به‌دست آمدند) با داده‌های 
پیش‌بینی شده توسط نرم‌افزار از طریق آزمون ] انجام گرفت. نتایج 
این مقایسه برای مرحله گرم کردن و سرد کردن به‌ترتیب در شکل‌های 


٩و‏ ۱۰ ارائه شده است. مقایسه نتایج بدست آمده از مدل دارای سرفضا 
و مدل فاقد سرفضا نشان داد اختلاف معنی‌دار بین دو مدل طراحی 
شده در سطح احتمال یک درصد وجود ندارد. اما داده‌های دمایی حاصل 
از مدل دارای سرفضا مطابقت دقیق‌تری با داده‌های تجربی داشتند. 
مقدار میانگین ريشه مربعات خطا مربوط به مدل دارای سرفضا در 
مقایسه با داده‌های تحربی برابر با ۲/۲۱ بود. 


۹۶ نشربه پژوهشهای علوم و صنابع غذایی ابران. جلد ۰۲۰ شماره ۱. فروردین-اردیبهشت ۱۴۰۳ 


مدل دارای سرفضاً 
6 ۵۵0 ۳/1۱۲ ۷۲۵061( ۱۴۰ 
مدل تجربی 00061 )صعصنته۲ / و ها تا ۱۲۰ 
اتوکلاوا56101 / ۱۰۰ 
و 3 
ط 
2- 
8 
۰ ام 
3 
۰ مب 
1 ۲ ۰ 
۰ ۸ ۰ ۶ ۰ ۲ ۰ 
زمان 
(صتجص) عحطرزً ]" 
شکل ۷- نمودار دما در محل هندسی قرارگیری ترموکوپل در مدل تجربی و مدل شبیه‌سازی شده دارای سرفضا در عصاره مالت و نمودار 
دما در محیط اتوکلاو 
۲ 21 5۵266 6۵0 ۱6۵ (06م رم ممتاتومص ۲معصهی مط) اه 6ااههصرصع) اصمصتهمی مه ۳۳۵۵۵۲۵۵ 7۰ ۲۱۵۰ 
۰ 122 
< ۱۸ 
د ۱۶ ِ_ 
سس | ۲ از 
9 3 
۱ ض و 
م۲۵ | تجربی ۵ ۸ ۸ 3 
۱ 4 از ۶ 9 
40 | پیش بینی شده 4 2۵ ۴ 
ی ی تن 
سس ۱ 
1 ل ل ل ل ۱ 1 1 1 1 1 0 1 1 1 1 1 ] ‌ 


محر[ 


شکل ۲0-۸ در مدل پیش‌بینی شده و تجربی دارای سرفضا 
6 ۷۷۱۲۴ 001 مه )ممصنه مرن مه ۵016)۵0 ۵۲ عط) صا و 1۳ .8 بع۲1 


بهینه‌سازی و تصحیح ارزش پاستوریزاسیون 

با توجه به عدم وجود ۳0 استاندارد برای کنسرو عصاره مالت می‌توان 
مورد مقایسه قرار داد. در این پژوهش میکروار گانیس مکلستریدیوم بوتولیتوم 
به‌عنوان میکروب شاخص در نظر گرفته شد. با توجه به وجود هماهنگی 
تقریبی بین محتوای شیر تغلیظ شده و ترکیبات عصاره مالت به‌نظر می‌رسد 
0 ۳02۵ برای محصول مناسب باشد (1979 ,۲60820). با توجه به 
اينکه در فرآیند استریلیزاسیون از سیکل‌های حرارتی تعریف شده برای 
محصولات کارخانه چیکا استفاده شد میزان ۳0 محاسبه شده (۱۵/۲۷) 


ونر -1 


بیش از میزان مورد نیاز برای عصاره مالت بود و با توجه به اینکه میزان 
0 اعمال شده سبب حرارت‌دهی بیش از حد محصول می‌شود بنابراین 
می‌توان از سیکل‌های حرارتی ملایم‌تر استفاده کرد و زمان حرارت‌دهی را 
کاهش داد. به‌منظور بهینه‌سازی میزان حرارت‌دهی عصاره مالت» پروفیل 
حرارتی آن, اصلاح شد. روش گرافیکی سیمپسون" نشان داد با کاهش ۱۱ 
دقیقه از زمان فرآیند در انتهای مرحله حرارت‌دهی می‌توان به نتیجه 
مطلوب (۳۵۵/۲۵ ) رسید ( ,.اه ۶ 51050 :2004 .۵1 6۶ 5150 
شا 19 تموتار کیرات عمایر محصول طن ونان راتفر مال 
تجربی و مدل اصلاح شده نشان می‌دهد. 
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دمای آزمایشگاهی 


۵ همم مه م)ر۲ 
شکل -٩‏ دمای پیش‌بینی شده به دمای آزمایشگاهی در مختصات (۲/۷۰- و۰ و۰) در مدل دارای سرفضا از آغاز تا پایان مرحله حرارت 
دهی 
۶ 00 6ظ) [ذاصت 5۵26۵ 620 ۱۷16۵ ]۵۵06 1 (2.76- و0 و0) 26 6اه صرع) همه مه و۱ اه مره ۳۸۲۵۵0۵۱۵۵۵۵ 10۰ .۲12 
۱2۵( 


۷ < 1.0155: - 06 
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دمای آزمایشگاهی 
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شکل ۱۰- دمای پیش‌بینی ده به دمای آزمایشگاهی در مختصات (۲/۷۲۱- و۰ و۰) در مدل دارای سرفضا در مرحله خنک کردن 
عصنا00 6ظ) رز مدصو 620 16 ۵001 (2۰76- و0 و0) 2 )همه همهم ۱ هه مره ۳۲۵۵0۵۲۵۲۵۵0 10۰ .۲12 


نتیجه گیری 

با توجه بهاینکه یکی از چالش‌های فراینداستریلیزاسیون کنسروها 
تولید محصولی سالم» یمن و در عین حال دارای کیفیت تغذیه‌ای 
تلوب اسهم تضییی سوت افقان خرازه یرای برمینن بت فرط 
بهینه فرایند و حفظ کیفیت محصول ضرورت دارد. همبستگی بالا بین 
دماهای پیش‌بینی شده و اندازه‌گیری شده. نشانگر کارآمد بودن مدل 
توسعه داده شده در پیش‌بینی دما در نقاط مختلف کنسرو عصاره مالت 
بود. همچنین با استفاده از پروفیل‌های دما در مدل» محل ناحیه کند 


۱ 


حرارتی تخمین زده شد. اطلاعات حاصل بیانگر وجود ناحیه کند حرارتی 
در یک سوم انتهایی کنسرو و نشانگر وجود جریان‌های جابجایی در 
محصول بود. با توجه به اينکه ناحیه کند حرارتی دیرتر از سایر نقاط 
گرم می‌شود لذا در صورت تأمین کشندگی لازم در اين ناحیه می‌توان 
از استریل شدن سایر نقاط اطمینان حاصل کرد. با توجه به میزان ۳0 
محصولات مشابه و استفاده از روش گرافیکی سیمپسون مقدار 0 
اصالاع تور ناه کند خرارتن کر فول کنسوی عضاره مات ۸/۲۵ 


محاسبه شد. این تحقیق برای بهینه‌سازی فرآیند استریلیزاسیون صنعتی 


۹۸ نشریه پژوهشهای علوم و صنایع غذایی ابران. جلد ۰۲۰ شماره ۱. فروردین-اردببهشت ۱۴۰۳ 


عصاره مالت با توجه به دما و زمان قابل استفاده است و می‌توان با توجه استریلیزاسیون را پیش‌بینی کرد. با بررسی سرعت انتقال حرارت می- 
به نوع غذای داخل بسته و ویژگی‌های آن زمان لازم برای توان موقعیت ناحیه کند حرارتی را پیش‌بینی و مدت زمان لازم برای 


استریلیزاسیون ناحیه کند حرارتی را محاسبه کرد. 
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شکل ۱۱- تغییرات دمای عصاره مالت طی زمان در مدل تجربی و مدل اصلاح ده 
هم 0:1160مه 0صد )هرن عمط ما مصت) مصتسا مرتاری اامص ۵ عمجم م۱6۵۵ :11 .م۲1 


1 0]6۳۵۵۵۵ 


204 ۱6۵۵0۲۵0۲۵ حتفم ۵۴ ممتهااصته آهعزمصصت ۲ 2002(۰) 2۸۰۲۸۰ بحعط »ی .۱۷۲۰۷۲ ب۲2۲1 ,.ت).ظ متصقدات هار 
۵66۳۵ ۳0۵ 0 09۳۱۱۵ .وعتصصهه ول متا آمطمناهانام‌صم مصتودا ماه 50۵112 عصتال صفی [)ط120عمظ 2 طا فع۳0111م 
005://001.018/10.1016/50260-8774)01(00039-5 .77-83 ,(1) 51 

اقصصمط) عمط 0۶ ممتاهعتاوع۷ 1۳ افاممحممورم 24 1۳6۵016۵06۵1 2002(۰) 2٩.۲2۰‏ یقعطن عک ,.۱۷۲,۲ ,نت۲2 ,.ت).ظ متصقطان اضر 
221-۰ ,(3) اد ,۱8۵66۱8 ۲۵۵0 ۵۴ ام .عمطمتاه 1۵00 ها دناتطام محصرعطان‌نهماه فسااتعدظ ۵ ممتاه۷ 1۳8۵001 
1016/50260-8774)01(00060-7 2/10۰ ۱۱۲05://001.0۲ 

۵ ۴ 06200۷21100 ۵۴ ممتاحوتاوع18۳۷ مه ,(1999) .۲ رولعقطامن ی .ابر بتعج ,۱۷۲۰۱۷ م۲2 ,.ت).ظ متصقطان اضر 
۵ 0 0۱۲۵ 0۳۲۱۲(۰) وعتصصمطرل مب امممتامانمهمی. عصلفی ممتاهمتاننماو عصتیت ‏ ۵00 8تبامنا 0مصصوع. 2 ظ1 
05://001.018/10.1016/50260-8774)99(00123-5ظ 207-214 42 روععاعو ۲ 

ما هه ۵ طمتامناصزی لهع1عممصیی ۲ ,(1999) .ظ وولتهطمن عک ر.ما.ع2 معط ,]۷۲۰۷ ,۲2۲1 ,.ت).ظ متصعطام ار 
,55-4 ول 8۵۱۱۵۵8۵ ۳۵۵0 ۵۴ سول .ععتصهدر . لت تقصمتاماه مهم وه 1۵00 0عصصق ‏ ۵۶ . عصتاهعط 
۱۲05://001۰0۲2/10۰1016/50260-8774)89(000735-4 

2013(۰) ,۷۲ مصتَلع۷11 عک بب م1125 ۷ ویژ.ل مهتصتا0 پیکل.]۷ مقلفتهت ر.ط ,0طهع-182ظ ویظ ,وع۷680-عع1تظ ر,.ضیظ ر0عطماظ 
,(32)2 ,0۳1۲0 ۳۲۵۵۵ 1۵0۵05۰ ۵۵2012826060 ط0ا۱۵۲۱۱122٩‏ لمصصعطا مه ۵۶ آمعجمی ادن ۶0۲ ممتادمنصتامون محصتا لقم۴ 
82/10/2 05://001.07ظ :392-403 

٩016006۰ ۰‏ 0و۸ .0۵6۲۵1۲۵7۵ ۵18۵۵۳8۵ ۲۵۵۵ -(1979) .۸۵ ,صحصصعءظ 

عص هط حعصه 0عمصنصظ طم2م گم ومزعماونط مساههم‌ما ۵ ممتاهایصته لیمیا 2021(۰) .۲ رماع کی رب یله 
۰ ,63 رععع6۵ 1۱۵۱۵ ۵ هل ۲۱۵۱۱۵۵09۵ .ماسقا متصصاه . همست فاصعتومه. نق. اصمع 
68 1 6/1 10.101 /05۹://001.0۲2) 

۲ 0660 عصتییل فععصمط انلمتان 4ص عصتاه0مصه تعتفصه) وققجه هه ۱۵۵۲ .(2022) .۸۵ بلته/۲۱۳ عک ,بآ بفحطتقطه ریک رطوجر1 
0۵:00 43)4(۰ ,هو ووعع۳۵ظ 0۴1۳۵۵۵ 0۵۱۵۱۵ ۲۵۷16۷۷۰ 51۷۵صعطه ۵۲ ۸ ۲۳۱۷18۵۰ 
فصه 10۱2/1 1 عع1وصهتا عمط ۵ فامعه ممتاع)۲0 220121 ۲۵۲ رای تمطمتامادام‌صصمی خر .(2012) ۱۷۲۰ وتهاتا 1 عک ,.۲ بتاع۲۲008 
,1۲۵۱۹/۵1۰ کو ۷۵ ۵۵ 01 ]۵ ۵۱۵ ۲۱۲۵۱۱۵۵0۳۸۵ 82050۵66۰ هه صصعاونره 1000 نهد تمه و11 عصتصتهاوهم 
۲20916۲.2001.1) دوعص معط ۰1 [/008://001.012/10.1016 .3774-3788 ردد 

کته رتم۷2 متصان ۵ ناه 0۵0۴۵۵۵۲۵۵۷۵ 2018(۰) ۲۰ ۱۷۲۱۲9۹۵۵۵02 ی بط رتطاک۷1 127 ر.9 بتل۷0056 بیط ,۲۵02۷1 
561۵166 ۳۵۵ ۲۱۱۱۵۷۵۷۵ .608۵60۵6 مفقصطها ۵۶ بزاتلدتن فطع من فمطامصظ مصتاجعط نامه لمممتامه همم 220 
4 /۸05://001.078/10.1016 .225-2300 و47 رععآع661:110108 1 ۲۱۱۵۲8۱۳۱8 

مصصتا جمتاجمتاتهاه ۵۶ جمتاقصتایی معط ۵۲ میامتصطهعا لمطمتامانم‌جم «ع خر 2004(۰) .۸.0۲ رتصجدام) اتلطاه ک ,۱۷۲۰۷۲ ,۲۵۲10 
7 ۳))0۹://001,0۲8/10۰۱016/.060.2003 ۰ 323-53 ,43 رهااوععع ۵( 0۳00 وهع و1 0۱۵6۵1 .1000 08عصصهه ظ1 
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۶ عاتاتاوص1 .ع0مطاعصط ]وم صرح فصمتاهم0عم؟ ]مها ]۱21 مات ,3897 ۱۱۲۷۲8۱۲ 191۳۴ جوا ۵۶ متاطانامرف۴ متصصهاو[ 
۷۰ 0۳۴ ۳۲1,۲6 ۴۳۳۲۲ ۵۸۱۷۲۲۵ ]19 .صف ۵۶ طمنهه‌ومک امتتاونصا مه وهعهلصها 

۱ ۸۱۱۵/۱۱۵۵ .9م0۳2 برطه‌طانقی ۵۶ طمتاحصتصعاع0 10۲ 0مط)عصظر محرمتطصنح طا خصفجصع۷متمرحصاً :(1952) ,۲.۸۵ ,1.0605 
۳)05://001.0۲2/106600612050 .219 ,(24)1 

معط عصتتل عمهمولهم وه 1۵۲ فاصمله‌تللهم بوصم کقعط 66001۷۵ جح گ۶ه ممتامصتصهاه(1۲ 2007(۰) 1:0۰ ,عصصهطم]۱۷ 
5-.16/[.10006۳82.6 10.10 /01005://001.0182 1166-1171۰ ,79 ,۱۱۵۱66۲۸8 0۳۵۵۵ ۱۱۵1 ول .عوعع۲۵ ٩611112۵109‏ 
۰ ( ۵۵ ما مصته ۵ عصتاملمصه تعامصه اهمز۲ .(2021) ۲۰ بتصطه009۳ عک ر.ظ متکل0۳00 .۷ بتعظ میک ,۱2۲۹2 
44)11(۰ ۵166۳8۵ . ۳۵6۵55 . ۳۵۵0 0۴ . 0۱1۱۵ م9۹ ص682 مان ...1۳01۷10۲۵1 
8 0005://001.072/10 

فصملوز1] ,۷۲ ,۱۷۲ ,بآ :1 60-81 ۳۳ .و7۱۵0 ووع7۵6 ۲۵۵۵ 2001(۰) .۱۱ رکامط0۵۵۲ عک ر.ظ بمه۷۵۲۵0۷ ,.]1۷.ظ ۱16012 
۰ 2۵ .۲۲688 تک .۲ع]قصه: 1 ووه]۱۷ 0ص )۳۵۵ ۶ه عصذاع۱۷]۵0 16 ,۳۱۵۲۲۵۵ ,1۷۲ .۲ .۸۵ .]۷۲ 200 

مج جونوع0 )صمتته 2۵0 6116001۷6 مج - (بالن) ومتصمورل مت لمممتاهانم‌جمن 2006(۰) 1.۷۷۰ مصتا کی با صماتول 
600-۰ ,(7)11 رآ 1661010/0 عک ۵62066 ۳۵0۵0 7 ۲۳۵/۵5 ۲۵۷۱۵۷۰ ۵ :1۳0090۳۷۰ 1000 عطا 10۲ 01م1 و791 202 
۰.4 10۰1016/12 /۱005://001.0۲2 

۶ مهطهمه]۱۷ ۱۱۵۸۵۲۱۵۱۰ ۵۵0 ک مایت وه ۵ و20۵۵ و02 2012(۰) ۰ رتتج‌طاله ‏ ره 1۷۲ ,۲۹۵۸2۵۷1 
۰ ,۲1۷۵۲۹16۷ 

۲۲۵۲۷۵۵۱ 160۵1۵21 )فع۴ ۱۷۱۵016 .و۲۵۵0 ۵۴ جع۴۳۵۵۵۲۲ ۳۳۹۱6۵ 2006(۰) .9.0 راصحنا عک رو مطنطه 
۷۵۱2 همه مد بامممصم ماع ام0ممط لهمتاممصمط ۱۷2۵ 2004(۰) ۱۷۲۰ بامطام۷ عک ر.ک رتفقطمصله ویک بطم0وموصل 
٩۳3206. 09۲۱۱۵ 0 ۵0‏ حصباوداظ ععمع صاً ۵000و ط 0ععامده دمط‌امص قاصهاعماه؟ :صمتام2تهتامه مصلودع0۳06 2۳20 
05://001.018/10.1016/50260-774)03(0024-72 :153-162 ,63 رو۱۱6ع۳۲۱۱ 

صوع0۳0 ۲۵۱0۲ مه ممتاه‌اناصصته تعتفصهط ۲۱۵2۶ 2015(۰) ,۱۷۲ رمهل‌صفط(مطعظ تعطمطه عک ولا رتهصتصه2 بلط مطفطلصماه ۱۷ 
۰ ا ص21 1۲1210 تمد ا مرمعآه۲۱ ۵2860 7/۵۵۵۲ ۵۶ عطلووعع0۳0 لقصصعطا عم اممصامن ره 


۹۹ 


12 


13. 


14 


15۰. 


16۰ 


17 


18۰ 


19۰ 


20. 


21۰. 


